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We present measurements of the J/tp invariant yields in ^/s NJV =39 and 62.4 GeV Au+Au collisions 
at forward rapidity (1.2 < \y\ < 2.2). Invariant yields are presented as a function of both collision 
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centrality and transverse momentum. Nuclear modifications are obtained for central relative to 
peripheral Au+Au collisions (Rcp) and for various centrality selections in Au+Au relative to scaled 
p+p cross sections obtained from other measurements (Raa)- The observed suppression patterns at 
39 and 62.4 GeV are quite similar to those previously measured at 200 GeV. This similar suppression 
presents a challenge to theoretical models that contain various competing mechanisms with different 
energy dependencies, some of which cause suppression and others enhancement. 



Heavy quarkonia are bound states of charm-anticharm 
or bottom-antibottom quarks. It was proposed over 25 
years ago that these states would be color screened in a 
quark-gluon plasma (QGP), thus suppressing their final 
yields in relativistic heavy ion collisions (|. The NA50 
experiment at the CERN-SPS measured a significant sup- 
pression of J/ip and ip' in Pb+Pb collisions at v / s^" = 
17.2 GeV, which was interpreted as indicating the onset 
of quark-gluon plasma formation [|J. However, measure- 
ments by the PHENIX experiment at the Relativistic 
Heavy Ion Collider (RHIC) indicated a similar level of 
nuclear suppression at midrapidity in Au+Au collisions 
at = 200 GeV @- Additional PHENIX results 

also indicated a larger suppression at forward rapidity 
1.2 < \y\ < 2.2 compared with midrapidity, despite the 
expectation of a higher energy density and temperature 
for the medium at midrapidity. Perhaps more surprising 
is the comparison of the recent higher statistics PHENIX 
forward rapidity J / ip suppression g and the ALICE ex- 
periment measurement in Pb+Pb at 2.76 TeV || at the 
Large Hadron Collider (LHC). These results indicate sig- 
nificantly less suppression for the most central Pb+Pb 
events at the LHC compared with Au+Au events at 
RHIC. Results at RHIC and the LHC at larger transverse 
momentum (pr > 4 GeV/c) suggest the opposite, 

with more suppression at the LHC compared to RHIC. 

These measurements contradict an interpretation 
based solely on color screening, and require the influence 
of other physics. There is an additional class of effects re- 
ferred to as "cold nuclear matter" (CNM) effects that are 
not due to the creation of a hot medium and thus can be 
probed via p(d)+ A collisions. These CNM effects include 
the modification of the initial incoming flux of quarks and 
gluons in the nucleus as described by nuclear-modified 
parton distribution functions (nPDFs) ||, breakup of 
the J/tp precursor cc state while traversing the nucleus, 
and initial-state parton energy loss. CNM effects have 
been studied in detail in p+A collisions at <Js NN = 17-42 
GeV pHBI, in d+Au collisions at ^/s^/ = 200 GeV by 



the PHENIX experiment |lj-|19(l , and p+A results from 
the LHC are anxiously awaited. In addition, there may 
be effects in the QGP other than color screening. These 
include the possible coalescence of originally uncorrelated 
c and c quarks or the recombination of breakup c and c 
pairs resulting in a competing enhancement effect (see 
for example Refs. [EOl Ell ) . This coalescence effect is ex- 
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pected to grow as the density of c and c increases. A 
recent review of many of these phenomena is given in 
Ref. [||. 

All of this highlights the importance of measuring J /ip 
and other excited quarkonia states over a broad range 
in ^Js NN \ thus varying not only the temperature of the 
medium, but also the c and c production and the cold nu- 
clear matter effects. In this paper, the PHENIX collabo- 
ration presents first measurements of invariant yields and 
suppression for J/ip at forward rapidity 1.2 < |y| < 2.2 
in Au+Au collisions at y/s NN = 39 and 62.4 GeV. 



I. DATA ANALYSIS 

The PHENIX experiment collected data in 2010 for 
Au+Au collisions at y/s^ = 39 and 62.4 GeV as part of 
the RHIC Beam Energy Scan program. After good run 
selection cuts, the data set includes 2.0 x 10 8 events at 
39 GeV and 5.5 x 10 8 events at 62.4 GeV. The PHENIX 
experiment is described in detail in Ref. [^3|. The J/ip 
measurement at forward rapidity is made via the dimuon 
decay channel with two forward angle muon spectrome- 
ters, as detailed in Ref. The muon spectrometers 
have acceptance over the range 1.2 < \<q\ < 2.2 and over 
the full azimuth. The two spectrometers comprise an ini- 
tial hadronic absorber followed by three sets of cathode 
strip chambers which are inside a magnetic field, referred 
to as the Muon Tracker (MuTr), and then five planes 
of Iarocci tubes interleaved with steel absorber plates, 
referred to as the Muon Identifier (MuID). Muon can- 
didates are found by reconstructing tracks through the 
magnetic field in the MuTr and matching them to MuID 
tracks that penetrate through to the last MuID plane. 

In Au+Au collisions at y/s~^ — 39 and 62.4 GeV, the 
events are selected with a minimum bias (MB) trigger 
utilizing the Beam-Beam Counter (BBC). The BBC com- 
prises two arrays of 64 quartz Cerenkov counters covering 
pseudorapidity 3.0 < \rj\ < 3.9. The MB trigger requires 
at least two hits in each of the BBC arrays and a recon- 
structed collision z-vertex of \z\ < 30 cm, where z = is 
the center of the detector. The BBC total charge is used 
as a measure of the collision centrality (the impact pa- 
rameter of the Au+Au collision is monotonically related 
to the average total charged particle multiplicity). Fol- 
lowing the procedure used for Au+Au at ^fs^ = 200 
GeV, for each centrality selection the average number 
of nucleon participants ((Ar part )) and the average num- 
ber of binary collisions ((V co n))are estimated using a 
Glauber model of the collision p5f and a negative bino- 
mial parametrization of the charged particles per pair of 
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participating nucleons. The total fraction of the Au+Au 
inelastic cross section measured by the MB trigger is de- 
termined to be 85.7 ± 2.0% and 85.9 ± 2.0% at 39 and 
62.4 GeV, respectively. The minimum bias sample is di- 
vided into exclusive centrality bins that are categorized 
via the Glauber model comparison to the BBC charge 
distribution as given in Table |. Note that the centrality 
selections used here are wider than in previous analyses 
for Au+Au collisions at y / s^7 = 200 GeV due to the 
smaller statistical sample of J / ips. 



TABLE I: Mean iVpart and JV co n values and systematic un- 
certainties in each centrality bin for Au+Au at 39 and 62.4 
GeV. 



y^GeV) Cent(%) 



(A/part) 



(iVcoll) 



30 



62.4 



0-40 
40-86 

0-20 
20-40 
40-60 
60-86 



204.4 ± 4.4 
34.1 ± 1.6 

274.8 ± 3.8 
138.7 ± 4.7 
59.7 ± 3.9 
14.7 ± 1.2 



444.8 ± 50.3 
43.5 ± 3.7 

689.9 ± 78.9 
270.5 ± 27.5 

85.7 ± 9.1 
14.3 ± 1.7 



For each centrality selection and beam energy, we ex- 
tract the number of J/ip counts following a method iden- 
tical to that used in Ref. Q. All unlike-sign muon can- 
didates are paired to calculate the invariant mass distri- 
bution. Underneath the J/ip signal are continuum back- 
ground counts both from uncorrelated tracks and from 
correlated physical backgrounds such as open charm de- 
cay (e.g. DD where both decay semi-leptonically to 
muons), open bottom decay, and the Drell-Yan pro- 
cess. First the uncorrelated background is calculated 
via an event mixing method with pairs from different 
Au+Au events with the same centrality and z-vertex. 
This background is then normalized using a compari- 
son of real-event and mixed-event like-sign pairs. Af- 
ter subtraction of the uncorrelated background, we fit 
to the remaining correlated dimuon spectrum with an 
acceptance-modulated J/ip line shape (determined from 
a full GEANT |26| simulation of the PHENIX detector) 
and an exponential folded with the acceptance to model 
the remaining correlated physics background. Utilizing 
different assumptions about the line shape, different un- 
correlated background normalizations, and different in- 
variant mass ranges for the fit (as detailed in Ref. [Q), 
we determine the systematic uncertainty on the extracted 
J/ip signal counts. The total J/ip sample corresponding 
to all centralities is approximately 170 counts at y/s NN — 
39 GeV and approximately 1060 counts at ^s NN = 62.4 
GeV. The invariant mass distribution of unlike-sign pairs, 
mixed-event pairs, and the subtracted distributions are 
shown in Figure [l]. The signal extraction procedure is 
quite robust and the systematic uncertainty is of order 
2-10%. 

The J/ip invariant yield is expressed as: 




mass (GeV/cf) 



FIG. 1: (color online) The unlike sign invariant mass distri- 
bution (blue) for all Au+Au centralities and both muon arms 
is shown along with the uncorrelated background calculation 
from mixed event pairs (red) for 62 GeV (top two panels) and 
39 GeV (bottom two panels). The lower panel in each pair 
shows the subtracted distribution and a fit to the data. 
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where _B MM is the branching fraction of J/ip to muons, 



V 



J/iP 



is the number of measured J /ips, Nevt is the 



number of events in the relevant Au+Au centrality se- 
lection, Ae is the detector geometric acceptance times 
efficiency, and Apr and Ay are the bin width in px and 
y, respectively. For the pr-fntegrated bins, we similarly 
calculate B^^dN/dy = Nj/^/(AeNEVT&y)- We evaluate 
the acceptance and reconstruction efficiency by running 
PYTHlA-generated |2?j] J /ips through the GEANT simula- 
tion of the PHENIX detector and then embedding these 
simulated hits into real Au+Au data events. These sim- 
ulated events are then reconstructed using identical code 
to that used in the real data analysis, and the overall 
acceptance and efficiency (Ae) is determined for each 
Au+Au centrality selection. In measurements at higher 
energy where the multiplicity is larger, there are large 
drops in the efficiency for more central collisions; but for 
these lower energies, with their lower multiplicity, there 
is no significant loss of efficiency for central collisions. 
There is an additional check on the efficiency of each 
MuTr and MuID plane that is determined via a data- 
driven method. The invariant yields are calculated sepa- 
rately for each of the two muon spectrometers and then 
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a weighted average taken. These results agree within un- 
certainties in all cases. 

Two categories of systematic uncertainties on the in- 
variant yields are shown in Table |n]: type A are point- 
to-point uncorrelated, and type B are correlated (or anti- 
correlated) point-to-point. The uncertainties listed in or- 
der are from uncertainties on the J /ip extracted yield as 
described above, the detector acceptance, the acceptance 
and efficiency over the rapidity range f.2 < \y\ < 2.2 
from the assumed PYTHIA input distribution, the abso- 
lute check on the MuTr and MuID hit efficiencies, and 
the matching of dead areas in the real data and GEANT 
Monte Carlo (MC) simulation. 



TABLE II: Systematic uncertainties 



Description 


Contribution 


Type 


Yield extraction 


2-10% 


A 


Detector acceptance 


5% 


B 


Input y,pr distribution 


4% 


B 


MuTr efficiency 


2% 


B 


MuID efficiency 


4% 


B 


DATA and MC mismatch 


4% 


B 
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FIG. 2: (color online) J/ip invariant yields (scaled by 
l/(^Vcoii}) are shown for Au+Au collisions at 39, 62.4, and 200 
GeV as a function of the number of participating nucleons. 
The solid error bars represent the uncorrelated point-to-point 
uncertainties (quadrature sum of statistical and type A); and 
the boxes represent the correlated (type B) systematic uncer- 
tainties. 



II. RESULTS 

Figure || shows the final calculated J/ ip invariant yield 
integrated over all px in y / s^" = 39 and 62.4 GeV 
Au+Au collisions as a function of centrality, categorized 
by the average number of participants (N paT t). The yields 
have been rescaled by 1/ (N co \\) . For comparison, the pre- 
viously published J/ip invariant yields in the same rapid- 
ity range 1.2 < \y\ < 2.2 from ^/J^ = 200 GeV Au+Au 
collisions are also shown j|. The vertical error bars are 
the quadrature sum of the statistical and type A sys- 
tematic uncertainties, and the boxes represent the type 
B uncertainties. As expected, the J/ip yield is larger in 
Au+Au collisions at larger center-of-mass energy. In ad- 
dition, the yield per binary collision is decreasing with 
(N paT t) at all three energies, indicating increasing nu- 
clear suppression for more central collisions. Figure ^ 
shows the invariant yield as a function of px , plotted at 
the center of each px bin, for y/s NN — 39 and 62.4 GeV 
Au+Au collisions. 

The nuclear modification of J/ip yields can be cate- 
gorized in different ways. Because the PHENIX experi- 
ment has not yet measured the p+p reference baselines 
at y/s NN = 39 or 62.4 GeV, we first discuss the J/ip Rcp, 
the nuclear modification of central relative to peripheral 
classes of events as defined below: 



Rcp = 



dN t; Ajdy (central) 
dN/ ^i Au /dy (peripheral) 

W *coll) 



(2) 



10 = 



a io-< 



10- 



10' 
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FIG. 3: (color online) J/ip invariant yields in minimum bias 
Au+Au collisions at 39 and 62.4 GeV as a function of trans- 
verse momentum. The solid error bars are the quadrature 
sum of the statistical and type A systematic uncertainties, 
and the boxes represent the correlated (type B) systematic 
uncertainties. 



The Rqp values are shown in Figure |J and in Table |II 
for Au+Au at 62.4 GeV. Note that the peripheral bin 
selection for Au+Au at 62.4 GeV is 60-86% centrality 
with a corresponding (N co \\) — 14.3 ± 1.7. Many uncer- 
tainties in the invariant yields cancel for Rqp and the 
dominant uncertainties are from the normalization with 
respect to the peripheral bin including the uncertainties 
in the (iV co rj) values for each bin. There is an additional 
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type C global systematic from the uncertainty in the pe- 
ripheral (-/V co n) value listed in the figure legend and in Ta- 
ble [Tn| the other systematic uncertainties are included in 
the boxes on each data point. For comparison, we show 
the published Au+Au results at 200 GeV || where the 
peripheral selection is 60-93%, with a quite comparable 
(Acoii) = 14.5±2.7. Within uncertainties, the centrality- 
dependent nuclear modification from peripheral to cen- 
tral collisions at the two energies are the same. 

TABLE III: PHENIX 39 and 62.4 GeV J/tp R C p vs Cen- 
trality with statistical uncertainties and Type A, B, and C 
systematics. 

y/s(GeV) Cent(%) Rcp Stat Type A Type B Type C 
39 0-40 0.554 0.112 0.028 0.138 0.047 



1.4 

1.2 : 
1- 



o°- 8 r 



0.6 

0.4 1 
0.2 : 



200 GeV, PRC 84, 054912 (2011) 

global sys. =±19.6% 

Peripheral (60-93%) : <N >= 14.5 ± 2.7 

62.4 GeV 

global sys. =±11.6% 

Peripheral (60-86%) : <N >= 14.3± 1.7 




J/y -» |in, 1.2 <|y| < 2.2 



50 



100 150 200 250 300 350., 400 

N part 



62.4 



0-20 
20-40 
40-60 



0.266 0.050 
0.353 0.064 
0.471 0.089 



0.005 
0.008 
0.013 



0.036 0.031 
0.045 0.041 
0.060 0.055 



For the Au+Au results at 39 GeV, the statistics do 
not allow any centrality dependence of Rcp and only a 
single value is calculated for the ratio between 0-40% to 
40-86% centralities, as shown in Figure | and in Table |n|. 
The published Au+Au results at 200 GeV are rebinncd 
to have a peripheral centrality selection of 40-93% to ap- 
proximately match the number of binary collisions for the 
peripheral denominator. Within uncertainties the results 
agree; however, the limited statistics in the Au+Au at 39 
GeV preclude any strong conclusions. 

The centrality dependence as quantified via Rqp is 
not a replacement for the nuclear modification factor 
Raa (relative to the p+p baseline) since J/ip yields may 
change already in peripheral Au+Au collisions, in partic- 
ular from cold nuclear matter effects. In addition, Rqp 
has significant uncertainties from the more limited statis- 
tics and the larger systematic uncertainty on (N co n) for 
the peripheral bin. The PHENIX experiment has no data 
for p+p collisions at 39 GeV, and only a limited data set 
was recorded during 2006 for p+p collisions at 62.4 GeV. 
However, p+p measurements do exist from fixed target 
p+A experiments near 39 GeV and from ISR collider ex- 
periments at 62 GeV. In the Appendix, we discuss in 
detail these results and compare them with theoretical 
calculations within the Color Evaporation Model (CEM) 
from R. Vogt p8 



29\ to determine a p+p reference. 
We quantify the nuclear modification factor Raa with 
respect to the p+p reference as follows: 



Raa 



1 dN AA /dy 
(T AA ) daPP/dy 



(3) 



where dN AA /dy is the invariant yield in Au+Au col- 
lisions, da pp /dy is the p+p cross section, and (Taa) 
is the nuclear overlap function (where (Taa) = 
(A r coii)/crj^ astic ). Unlike 200 GeV, the 39 and 62 GeV 



FIG. 4: (color online) J/tp R C p for 0-20%, 20-40%, and 
40-60% (central) relative to 60-86% (peripheral) Au+Au col- 
lisions at 62.4 GeV. For comparison, Rcp results from Au+Au 
collisions at 200 GeV are shown with a peripheral bin of 60- 
93%, where the {N co n) value is a close match. The solid error 
bars are the quadrature sum of the statistical and type A sys- 
tematic uncertainties, and the boxes represent the correlated 
(type B) systematic uncertainties. 



p+p references are determined from other measurements 
rather than being from our own, and systematic uncer- 
tainties will not cancel in the ratio. Our estimates for the 
J/4> P+P cross sections in the range 1.2 < \y\ < 2.2 for 
39 and 62.4 GeV are shown in Table |v|, and are detailed 
in the Appendix. The J/ ip Raa for Au+Au collisions at 
39 and 62.4 GeV is tabulated in Table |y| and shown in 
Figure |as a function of the number of participating nu- 
cleons (-/Vp ar t), along with the previously published 200 
GeV results Q. The type C global scale uncertainties, 
from the p+p references, are listed separately in the leg- 
end. At both 39 and 62.4 GeV, there is slightly less J/tp 
suppression than observed in Au+Au at 200 GeV. How- 
ever, particularly for 62.4 GeV, since we have no reliable 
p+p reference from our own measurements, the Raa re- 
sult could shift down by the quoted 29% systematic un- 
certainty, bringing the data into agreement with the 200 
GeV result. 



TABLE IV: Estimates used for the 39 and 62.4 GeV J/tp 
p+p cross sections along with their uncertainties. See the 
Appendix for details. 



v^(GcV) 


da pp /dy 


39 


2.91 ± 19% nb 


62.4 


7.66 ± 29.4% nb 
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FIG. 5: (color online) J/tp Rcp for 0-40% (central) relative 
to 40-87% (peripheral) Au+Au collisions at 39 GeV. For com- 
parison, Rcp results from Au+Au collisions at 200 GeV are 
shown with a peripheral bin of 40-93%, where the (iVcoii) value 
is a close match. The solid error bars are the quadrature sum 
of the statistical and type A systematic uncertainties, and the 
boxes represent the correlated (type B) systematic uncertain- 
ties. 



TABLE V: PHENIX 39 and 62.4 GeV J/tp Raa vs Cen- 
trality with statistical uncertainties and Type A, B and C 
systematics. 





Cent(%) 


Raa 


Stat 


Type A 


TypeB 


Type C 


39 


0-40 


0.439 


0.043 


0.020 


0.077 


0.083 




40-86 


0.793 


0.157 


0.011 


0.139 


0.151 


62.4 


0-20 


0.292 


0.039 


0.004 


0.042 


0.085 




20-40 


0.388 


0.047 


0.008 


0.056 


0.115 




40-60 


0.519 


0.067 


0.014 


0.073 


0.153 




60-86 


1.100 


0.150 


0.010 


0.155 


0.323 



III. DISCUSSION 

The collision energy dependence of the various compet- 
ing effects influencing the final J/tp yields are all quite 
different. Thus, the similarity of the J/tp nuclear modifi- 
cations Rc p and Raa from 39 to 200 GeV is a challenge 
for models incorporating the many effects. There was 
a prediction that the maximum J/tp suppression would 
occur near y/s^ — 50 GeV, as shown in Figure [[Jo| . 
As the collision energy increases the QGP temperature 
increases, and thus the J/tp color screening (labeled as 
direct J/tp suppression) becomes more significant. How- 
ever, in this calculation, the regeneration contribution 
increases with collision energy due to the increase in the 
total number of charm pairs produced and nearly com- 
pensates. This result is for J/tp at midrapidity and rela- 
tive to the total charm pair production (thus removing in 



FIG. 6: (color online) J/tp Raa at v /s^7 = 39, 62.4, and 200 
GeV. The solid error bars are the quadrature sum of the sta- 
tistical and type A systematic uncertainties, and the boxes 
represent the correlated (type B) systematic uncertainties. 
The global systematic uncertainties are quoted in the legend 
for each energy's results. 



this ratio possible changes in the charm pair production 
caused by initial state effects). 

Recently, the same authors have completed new cal- 
culations including cold nuclear matter effects, regener- 
ation, and QGP suppression specifically for J/tp at for- 
ward rapidity |3f], [32| . Figure |8| shows these results (in 
the so-called "strong binding scenario"). The contribu- 
tions of direct J / tp and regeneration are shown separately 
(and scaled down by x0.5 for visual clarity). The in- 
clusion of cold nuclear matter effects and the forward- 
rapidity kinematics slightly reverse the trend seen in Fig- 
ure and now the total J/tp Raa follows the ordering 
.Raa (200 GeV) < i?AA(62 GeV) < i? A A(39 GeV) (though 
by a very modest amount). Also shown in Figure |§| are 
the PHENIX experimental measurements that, within 
the global systematic uncertainties, are consistent with 
the theoretical calculations. 

These results highlight the need for p+p reference data 
at both 39 and 62.4 GeV from the same experiment. In 
addition, the cold nuclear matter effects are likely to be 
different at the different energies (an important input for 
the calculations in Ref. |5l]]). The x distribution of glu- 
ons for producing J/tp at 1.2 < \y\ < 2.2 changes as the 
colliding energy decreases. In a simple PYTHIA study, 
one finds that the average gluon X\ and xi for producing 
J/tp between 1.2 < \y\ < 2.2 is 0.14 and 0.01 for 
= 200 GeV, 0.32 and 0.03 for = 62.4 GeV, and 

0.43 and 0.05 for y/s NN = 39 GeV. The large uncertain- 
ties in the gluon nPDF for the anti-shadowing and EMC 
regions || leads to an additional ±30% uncertainty in 
the J/tp initial production for the central Au+Au case. 
Future measurements in p(d)+A collisions at these en- 
ergies are clearly required in order to reduce this large 
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IV. SUMMARY 
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The PHENIX experiment has measured the invariant 
yield of J ' at forward rapidity in Au+Au collisions at 
39 and 62.4 GeV. The nuclear modification, when formu- 
lated as i?cp (the ratio between central and peripheral 
event classes), indicates a similar suppression pattern at 
the two lower energies to that previously published for 
Au+Au collisions at 200 GeV. Using a p+p reference from 
other experiments and from a Color Evaporation Model 
calculation, results in an Raa with slightly less suppres- 
sion at these lower energies. These results are consistent 
with theoretical calculations dominated by the balancing 
effects of more QGP suppression as well as more J/ip re- 
generation for high-energy collisions. However, any firm 
conclusion regarding the overall level of suppression from 
the QGP requires additional p+p and p(d)+A data at 
these energies. 



FIG. 7: (color online) The number of J/ip per produced 
charm pair (xlO ) in Au+Au central collisions (Apart = 
360) at midrapidity. Shown are the direct J ftp and regenera- 
tion contributions. Calculation details and figure from fed]. 
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FIG. 8: (color online) The J/tp nuclear modification factor 
Raa as a function of the number of participating nucleons 
Apart for ys^T = 39, 62.4, 200 GeV Au+Au collisions. Cal- 
culation results are shown from j3lj for the total J ftp Raa 
and the separate contribution of direct J/tp suppression and 
regeneration (scaled down by x0.5 for visual clarity). The 
PHENIX experimental data points are shown for comparison. 



uncertainty contribution. 
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FIG. 9: (color online) J/ip cross section as a function of 
rapidity in p+p collisions at 200 GeV. The CEM calculation 
is shown as a black solid line with a gray band for its uncer- 
tainty. In comparison, PHENIX measurements are shown as 
red points 



FIG. 10: (color online) J /if) cross section as a function of 
rapidity in p+p collisions at 39 GeV. The CEM calculation is 
shown as a black solid line with a gray band for its uncertainty. 
Data points and fit function are the result of the p+A data 
extrapolation to p+p as described in the text. 



Appendix: Proton-Proton Reference 

In order to construct the p+p references at 39 and 62.4 
GeV, we utilize lower-energy data from Fermilab and the 
ISR, and also the Color Evaporation Model (CEM) cal- 
culations from R. Vogt ^] . These calculations have 
been extensively compared with J/tj) cross sections as a 
function of center-of-mass energy. First, shown in Fig- 
ure [9] is a comparison of the published PHENIX mea- 
surements for the J I ip cross section in p+p collisions at 
200 GeV H and the CEM calculation. For the CEM 
calculation, the solid line is the central value and the 
gray band represents the systematic uncertainty of the 
results. Using the same CEM framework, calculation re- 
sults for p+p at 39 and 62.4 GeV are shown in Figs. |l0| 
and 11 , respectively. It is notable that the predicted cross 
section at midrapidity drops by approximately a factor 
of 2.5 in going from 200 to 62.4 GeV, and then another 
factor of 1.9 in going from 62.4 to 39 GeV. The rapidity 
distribution also narrows as expected. 
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FIG. 11: (color online) J/ij) cross section as a function of 
rapidity in p+p collisions at 62.4 GeV. The CEM calculation is 
shown as a black solid line with a gray band for its uncertainty. 
Data points from ISR measurements are shown as detailed in 
the text. 



1. p+p at 39 GeV 

Fermilab fixed target experiment E789 ^5) has 
measured the invariant cross sections of J ftp in p+Bc, 
p+Cu, and p+Au collisions over a broad rapidity range 
at V s «« = 38.8 GeV. The rapidity coverage for p+Au 
was —0.1 < y < +0.7; and for p+Be and p+Cu was 
1.4 < y < 2.4. In addition, the nuclear dependence of 
the cross sections was measured by E866/NuSea jlEj and 
found to follow the functional form, o p +a = A a a p+p , 
where, 



a(x F ) = 0.960(1 - 0.0519^ - 0.338:4) (A.l) 

(as seen in Figures 2 and 3 of Ref. ) . 

Using this parametrization for the nuclear dependence, 
one can extrapolate versus A from the p+A J/t/j cross 
sections to those for p+p (A = 1) and obtain the p+p 
cross sections as a function of xp. After converting 
these to be versus rapidity, they are shown in Fig- 
ure [!(]. For the rapidity range 1.2 < y < 2.2 one ob- 
tains 2.91 ± 19%(syst) nb by integrating the fit func- 
tion. In comparison, the result from the CEM calculation 
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is 2.45^i'o 8 nb, which agrees well within uncertainties. 
Thus, we use this extraction from the experimental data 
for the 39 GeV p+p reference, as shown in Table IV. Sys- 



tematic uncertainties on this reference include 12% from 
the E789 p+A data and 15% to account for the quality 
of the fit and for its extrapolation in rapidity into the 
unmeasured 1.2 < y < 1.4 region. 



2. p+p at 62.4 GeV 



pidity coverage using the CEM calculation fitted to the 
Antreasyan measurement. For the uncertainties of this 
reference we take similar CEM guided integrals, but con- 
strained to the upper and to the lower limits of the that 
ISR measurement. This results in a 62 GeV p+p refer- 
ence of 7.66 ± 29.4% nb. We note that the midrapidity 
ISR points are somewhat low but nearly consistent with 
the CEM calculation, but since our data lies at large ra- 
pidity we rely on the Antreasyan ISR point. 



Experiments at the CERN Intersecting Storage Ring 
(ISR) measured the J/tp cross section in p+p collisions 
at 62 GeV |3(]-|38j. These results are shown in Ta- 
ble ^ and in comparison to the CEM calculation in Fig- 
ure [11 . Since our measurements lie in the rapidity range 
1.2 < \y\ < 2.2, the most important p+p measurement 
from the ISR for our purposes is that of Antreasyan jj6| , 
which covers a rapidity range of 0.89 < y < 1.82 and 
agrees quite well with the CEM calculation. Therefore 
we estimate the p+p reference by integrating over our ra- 



TABLE VI: ISR measurements of J/t/j in p+p collisions at 
62 GeV 



Reference 


Rapidity Range 


^j(nb) 


Antreasyan et al. J36| 


0.89 < y < 1.82 


9.21 ± 2.70 


Clark et al. §7| 


M<0.5 


10.2 ± 0.7 


Kourkoumelis et al. 


\y\ < 0.65 


14.8 ± 3.3 
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